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Abstract. Simultaneously with inventing the modern relativistic formalism of quantum electro- 
dynamics, Feynman presented also a first-quantized representation of QED in terms of worldline 
path integrals. Although this alternative formulation has been studied over the years by many au- 
thors, only during the last fifteen years it has acquired some popularity as a computational tool. I 
will shortly review here three very different techniques which have been developed during the last 
few years for the evaluation of worldline path integrals, namely (i) the "string-inspired formalism", 
based on the use of worldline Green functions, (ii) the numerical "worldline Monte Carlo formal- 
ism", and (iii) the semiclassical "worldline instanton" approach. 
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1. FEYNMAN'S WORLDLINE REPRESENTATION OF QED 

In 1950 Feynman presented, in an appendix to one of his groundbreaking papers on the 
modern, manifestly relativistic formalism of perturbative QED [1], also a first-quantized 
formulation of scalar QED, "for its own interest as an alternative to the formulation 
of second quantization". There he provides a simple rule for constructing the complete 
scalar QED S -matrix by representing virtual scalars and photons in terms of relativistic 
particle path integrals, and coupling them in all possible ways. Restricting ourselves, for 
simplicity, to the purely photonic part of the S -matrix (no external scalars), and moreover 
to the "quenched" contribution (only one virtual scalar), this "worldline representation" 
can be given most compactly in terms of the (quenched) effective action T[A]: 

TscalarfA] = / J 4 *Jzf scalar [A] = f ^ c~^ T [ ®x(t) e~ S ^ (1) 

J JO 1 Jx(T)=x(0) 

Here T denotes the proper-time of the scalar particle in the loop, m its mass, and 
Ix(t)=x(0) ^ X ( T ) a P at h integral over all closed loops in spacetime with fixed periodicity 
in the proper-time. The worldline action 5[jc(t)] has three parts, 

S = Sb + Sext + Sim (2) 

(see fig. [TJ. Of these, the kinetic term So describes the free propagation of the scalar, 
Sext its interaction with the external field, and Si nt the corrections due to internal photon 
exchanges in the loop. The connection to a standard Feynman diagrammatic description 
is made simply by expanding out the two interaction exponentials. 
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FIGURE 1. Perturbative expansion of the worldline path integral. 



The generalization of this representation to include multiple scalar loops and open 
scalar lines is straightforward. It yields a first-quantized representation of the full effec- 
tive action, and thus, by Fourier transformation, of the S-matrix. 

While for scalar QED this representation is essentially unique, when generalizing it 
to spinor QED one has a number of choices. First, worldline representations of spin 
half particles can be derived either from standard first-order Dirac theory, or from its 
second-order formulation (see [2] and refs. therein), based on the identity 



(0 + eA) 2 = -{d + ieA) 2 - l -ea^ v F^ v (3) 

Contrary to the situation in second-quantized field theory, in the worldline formalism 
the second-order approach is the more standard one, and I will restrict myself to it in 
this review (see EH) for the first-order approach). Using ([3]) one arrives at a worldline 
representation for the fermion QED effective action r sp i nor [A] which, at the one-loop 
level, differs from the scalar one (Q ) only by the addition of a global factor of — \, and 
the insertion of a spin factor S[x, AJunder the path integral 0, 



S[x,A] = tr r ^exp 



(4) 



Here trp denotes the Dirac trace and & path ordering. 

A more modern way of writing the same spin factor is in terms of an additional 
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S[x,A] = J exp — J dT^xif-xjf-iexi/^F^y^ 



(5) 



Here the path integration is over the space of anticommuting functions antiperiodic in 
proper-time, ^(^1)^(^2) = -^(^2)^(^1)' ¥^( T ) = -1/^(0). The main advantage 
of introducing this second path integral is that, as it turns out, there is a "worldline" 
supersymmetry between the coordinate function jc(t) and the spin function i/a(t) 0. 
Although this supersymmetry is broken by the different periodicity conditions for x and 
i/A, it still has a number of useful computational consequences. In particular, introducing 
a worldline superformalism allows one to combine the two path integrals, and to write 
down a formula for r sp i nor [A] which is completely analogous to ([I]), ([2]), including the 
internal photon corrections ifTTTl . 

A third, even more subtle, way of implementing spin on the worldline is the "Polyakov 
spin factor", a purely geometric quantity depending only on the worldline itself lfT2l[T3ll . 
This aspect of spin has been studied by a number of authors, but usually using the first- 
order formalism; the form of the spin factor appropriate for the second-order formalism 
has been established only recently lfT4l . 

A vast amount of work has been done on these QED worldline representations and 
their generalizations to other background fields and couplings (see lfT5ll for an exten- 
sive list of references). Nevertheless, most of the earlier work on this subject is con- 
cerned with formal aspects of relativistic particle Lagrangians, rather than with attempts 
at performing state-of-the-art calculations in quantum field theory (notable exceptions 
are US El and lfT8l0 . It is only during the last fifteen years that the usefulness of the 
first-quantized approach as an alternative to standard Feynman diagrammatic methods 
has been seriously investigated. This development was triggered by string theory, where 
first-quantized path integrals are the standard tool for perturbative calculations of scat- 
tering amplitudes, and by the fact that many amplitudes in quantum field theory can 
be represented as infinite string tension limits of the corresponding string amplitudes. 
Bern and Kosower lfT9ll investigated this limit in detail for the case of the QCD N gluon 
amplitudes, and found in this way a new set of rules for the construction of these ampli- 
tudes. Shortly later, Strassler's work ll20l showed that the same type of "string-inspired" 
representation of photon/gluon amplitudes can also be obtained more directly by evalu- 
ating the corresponding first-quantized path integrals using appropriate worldline Green 
functions. This approach was then generalized to some cases of multiloop amplitudes 
EHHDEZII23, as well as to QED amplitudes in a constant external field (MESl (see 
lfT51l for a review). Moreover, the success of this program led to the development of al- 
ternative techniques for the calculation of worldline path integrals. In this talk, I will 
restrict myself to the computational aspects of the worldline representation, and to the 
prototypical case, QED. I will discuss, in turn, three quite different methods which are 
now available for the computation of the scalar/spinor QED worldline path integrals, 
(i) the original "string-inspired" approach, (ii) the numerical "worldline Monte Carlo" 
method and (iii) the semiclassical "worldline instanton" technique. 



2. THE "STRING-INSPIRED" APPROACH 



The strategy in the "string-inspired approach" is simple. The path integral(s) will be 
manipulated into Gaussian form, after which they can be performed using worldline 
correlators with the appropriate periodicity properties. For the one-loop path integrals 
{j},© those are 1126112011 



</(Tl)/(T 2 )> = -G B (T h T 2 )S^ 
G B (T h T 2 ) 

(^(tO^(t 2 )) 

G F (T h T 2 ) 



PI - ^2 1 - — ( ^1 - ^2 



G F (z U T 2 )S^ 

sign(Ti - T 2 ) 



(6) 



Here (t) = x M (t) — x 1 ^ , where x 1 ^ denotes the loop center of mass, x^ = j J r dxx^ (t) . 
There is a certain freedom in choosing the coordinate correlator Gb (see, e.g., [15]); the 
one given above is usually the most convenient one, but some of the references given 
below use the alternative G#(ti, t 2 ) = \z\ — T 2 \ — (ti + z 2 ) + \i\i 2 . 
Now, to get gaussian form, 

• Expand all the interaction exponentials e~ s ^i x ( z )} etc. 

• For the effective action: Taylor expand the external field at the loop center of mass, 



A^(x(T)) = ^- d A^(x ) 
For the N photon amplitudes : Expand the field in N plane waves, 



(7) 



N 

AV(x(?)) = Y< £ i^ ki ' X{T) 

i=\ 



Exponentiate the denominator of the photon insertion terms, 



(8) 



87T 2 

df 



-Jo dXx L dx \ x ^)- x{ , 2 )f 



e 2 f°° dT f T f T 
~2 Jo (W/o dTl Jo ^i^O-i^exp 



{x{x x )-x{x 2 )f 
AT 



(9) 



For example, for the four-photon amplitude in scalar QED this procedure yields the 
following integral representation: 



Ql 
Ql 



Qf 



poo r j r r 4 pT 

r[*i,ei;...;*4,e4] = e 4 [4kT]- 2 e~ m T H dx { Q^ G ^ 

Jo T Jo 



Qa = Qt + Ql + Ql-Qf 



= Gb\2Gb23Gb34Gb4\Z4(1234) + 2 permutations 
= Gbi2Gb23Gb3iZ3(123)Gb4^4 'k + 3 perm. 

= G512G521^2(12) j 653^3 • k[(jB4j^4 ' kj + -653483 • £4 

+ 5 perm. 

= 65126521^2(12)65346543^2(34) +2 perm. 



(10) 
(11) 



653/^3 ' _ 654/^4 • /:/ 



(12) 



Z 2 (/j) = £/ • kj£j • fe; - £i • e/fc* • 

n r -, 

Zn{iii 2 ..A n ) = trf] ki j ®£ ij -£i j ®k ij (n>3) 
7=1 

(13) 

Here each of the integrals Jq d%[ represents one of the four photon legs moving around 
the scalar loop. Gsij stands for Ga(t/, T j) an d 65/y for its derivative. Repeated indices are 
to be summed over 1, . . . , N. In writing down the integrand, a number of integrations by 
parts have already been performed which eliminated all second derivatives of the G5//S 
that initially appear in the factor Q4. Moreover, this factor has been decomposed into 
terms Q™ according to "cycle content", indicated by the superscript, where a "cycle" 
is a factor of Gb^Gb^ Gb^ •••65.^. Such a 'r-cycle' always gets multiplied by a 
'Lorentz-cycle' Z n (i\i2- • .i n ), which is the trace of the products of the field strength 
tensors associated to the corresponding legs. The various Q™ 's are individually gauge 
invariant. This representation can be generalized to an arbitrary number of photons, and 
is unique if one requests manifest permutation symmetry in the photon legs [|T9ll20ll27ll . 
The generalization to higher loop photon amplitudes can be achieved either using ([9]), 
or by explicit sewing of pairs of external photons, or more directly using generalized 
worldline Green functions adapted to the multiloop graph topologies ll2Tll22ll28ll . 

In all cases the arising multiple integrals are equivalent to standard Feynman param- 
eter integrals. For example, in ( [T0| ) and its A/-point generalizations, the prefactor Qn 
relates to a Feynman numerator and the exponential factor to the universal denominator 
of one-loop A/-point integrals (see, e.g., Il29ll ). However, the string-inspired representa- 
tion has a number of interesting properties which are not usually manifest in standard 
Feynman parameter calculations: 



1 . The use of the cycle notation allows a more compact way of writing the N photon 
amplitudes than usual. 




FIGURE 2. Diagrams contributing to the three loop QED photon propagator. 



2. There is a simple "cycle replacement rule" which allows one to construct the 
integrand of the spinor loop N photon amplitude from the scalar loop one [|T9l 
120) . This implies that, in this formalism, the calculations of the same quantity in 
scalar and spinor QED are not independent; any spinor QED calculation yields the 
corresponding scalar QED result as a byproduct. 

3. The integral ( [T0| ) and its higher N generalizations represent the whole amplitude, 
with no need to sum over permutations. This property is not very relevant at the 
one-loop level, but becomes interesting at higher loop orders. In the QED case, it 
generally allows one to combine into one integral all contributions from Feynman 
diagrams which can be identified by letting photon legs slide along scalar/electron 
loops or lines. As an example, we show in fig. [2] the "quenched" contributions to 
the three-loop photon propagator. 

This property is particularly interesting in view of the fact that it is precisely this 
type of sums of diagrams which in QED generally leads to extensive cancellations 
between diagrams, and to final results which are substantially simpler than inter- 
mediate ones (see OO) and refs. therein). And for the two-loop QED j8 function 
indeed a way was found for calculating the corresponding integral in a way which 
avoided splitting up the multiple parameter integral into sectors with a fixed order- 
ing of the photon legs, and which led to dramatic simplifications ifTTTl . However, 
so far no generalization of the method used there to higher loop orders has been 
found. 

4. The string-inspired method provides a particularly convenient way of implement- 
ing constant external fields in QED calculations. Photon amplitudes or effective 
lagrangians in a constant field are obtained from the corresponding vacuum quanti- 
ties simply by substituting the vacuum Green functions Gb,f{^\ , ^2) by appropriate 
field-dependent Green functions ^b,f(^i^2' 9 F) 9 an d by a change of the free world- 
line path integral determinants [[24l [25]] (see also lUflh . In particular, the "cycle 
replacement rule" carries over to the constant field case. 

The efficiency of the technique for QED in a constant field has, at the one-loop level, 
been demonstrated by recalculations of the photon splitting amplitude in a magnetic 
field [32], and of the one-loop vacuum polarization in a general constant field 03) . At 
the two-loop level it has been extensively applied to the QED effective Lagrangian in a 
constant field. This includes recalculations of the standard two-loop Euler-Heisenberg 
Lagrangians (25l [3H [35), closed-form expressions for the weak field expansion coef- 
ficients of the magnetic two-loop Euler-Heisenberg Lagrangian |[36lh and a generaliza- 
tion of these Lagrangians to the case of a self-dual Euclidean field [37 j. I will show 
here only the last result, which is particularly nice: a Euclidean self-dual field fulfills 



Fjiv = ^e^vajS^ 7 "^' which implies that the square of the field strength tensor is propor- 
tional to the unit matrix, Fn V F v ^ = —f 2 $fi . Remarkably, this simplifies matters so much 
that all parameter integrals can be done in closed form, leading to the following explicit 
formulas for the two-loop scalar and spinor QED effective Lagrangians in a constant 
self-dual field 1571 . 

4inorW = -2«^^[3$ 2 (k) 

(14) 

2 

Here K = |^ and 

^^-S-^+s) <15) 

The simplicity of these results made it possible to use this self-dual case for studying 
the generic properties of the weak and strong field expansions of Euler-Heisenberg 
Lagrangians [38 j, as well as the corresponding maximally helicity violating components 
of the two-loop A/-photon amplitudes in the low energy limit Il37ll . 

Apart from the effective action in a constant field, the string-inspired method has also 
been extensively applied to the calculation of the full one-loop QED effective action in 
an arbitrary background field, using the derivative expansion [|39ll40ll4T1l . 

Moreover, at the one-loop level the method has been generalized to the finite temper- 
ature case, both to obtain representations of the N - photon amplitudes J42l |43) and for 
the effective action in a general background ll44ll . 

Finally, as should be clear from the above, work on the string-inspired technique in 
QED so far has been concerned mainly with purely photonic amplitudes. The extension 
to amplitudes involving also external scalars E3ll or fermions j45l|46l is possible without 
difficulties, although its practical usefulness is difficult to judge at present due to a lack 
of state-of-the-art applications. 



cfW ( K \ - a ^— — 



3. THE "WORLDLINE MONTE CARLO" APPROACH 

During the last few years it was found that Feynman's worldline representation ([T]) and 
its various generalizations are also very amenable to a direct numerical evaluation using 
standard Monte Carlo techniques 071 |48l |49]|. At the one-loop level, this approach 
has been shown to work well for QED effective actions in quite general background 
fields. This includes singular fields such as a magnetic field of the step-function type 
B71l . It also extends to the imaginary part of the effective action [50J, to be discussed 
in section 4 below. First steps towards a multiloop extension have been taken in ll5TTh 
although, as with all numerical approaches to quantum field theory, implementing the 
full renormalization program in such a formalism poses a formidable challenge. 



The Monte Carlo approach appears to hold particular promise for the calculation of 
Casimir energies for arbitrary geometries [|49l [52l [531 . Although this has so far been 
done only for scalar fields, not for the QED case, it shall be discussed here since it is 
interesting to see how Dirichlet boundary conditions can be implemented at the level 
of the path integral (see the recent Il54l0 for a treatment of boundary conditions in 
the "string-inspired" approach). For a (massless) scalar field in a background potential 
V(x), the gauge coupling term in Feynman's path integral ([I]) has to be replaced by 
— Jq drV(x(r)). Dirichlet boundary conditions on an (infinitely thin) surface Z can be 
implemented by choosing 



V(x) = X J^doS 3 (x-x G ) (16) 

with A — > oo. Considering the case of two disjoint surfaces E = Li U^2, it can then easily 
be shown H9l that the Casimir interaction energy between the surfaces is given by[] 



^Casimir = ~A= f ^ f d\ /^x(t) e^"* ©z[x(t)] 

2 V4^ Jo T2 J J 



(17) 



where @ £ = 1 if the path x(t) intersects both Ei and E2, and ©x = otherwise. The 
trivial time component of the path integral has been integrated out. xq denotes the loop 
center of mass of the remaining integral over spatial loops; the energy density at a point 
xo is obtained by restricting the path integral to worldloops with that point as their 
common center of mass. 

Figure [3] shows the Casimir energy density obtained by a numerical evaluation of ( fTTj ) 
for the case of two infinite perpendicular plates at separation a [52] 




FIGURE 3. Left panel: sketch of the perpendicular-plates configuration with (an artist's view of) a 
typical worldline that intersects both plates. Right panel: Density plot of the effective action density jSf 
for the perpendicular plates case; the plot shows ln(2(47la 2 ) 2 



When comparing with l49l note that they use a different normalization of the free path integral. 

Fig. 3 is reprinted from [ 52 ] with the permission of the authors and of the American Institute of Physics. 



A particularly nice example is the case of a sphere of radius 7? and an infinite plate at 
a distance a ll53ll . Here an exact result for the Casimir interaction energy is known as a 
function of a/R for a/R > 0.1 ll55ll . Fig. P] shows the result of the worldline evaluation 
for this case. As can be seen from fig. |5J the result is in very good agreement with ll55ll 
for the whole parameter range Fj 
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FIGURE 4. Contour plot of the negative Casimir interaction energy density for a sphere of radius R 
above an infinite plate; the sphere-plate separation has been chosen as a = R. The plot results from a 
pointwise evaluation of eq. flT] ) using worldlines with a common center of mass. 
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FIGURE 5. Casimir interaction energy of a sphere with radius R and an infinite plate vs. the curvature 
parameter a/R. The energy is normalized with respect to the zeroth-order proximity force approximation 
(PFA) [ 53 1 . The numerical worldline result is compared to the exact result of l55l and the PFA estimate. 



3 Figs. |4j and |5j are reprinted from [ 53 1 with the permission of the authors and of the American Physical 
Society (copyright by APS, http://link.aps.org/abstract/PRL/v96/e220401 ) . 



4. THE "WORLDLINE INSTANTON" APPROACH 



We come now to a third approach which is more specialized than the two previous ones, 
since it applies only to the imaginary part of amplitudes or the effective action. Although 
the basic idea was presented by Affleck et al. in 1982 lfT8lL it seems not to have been 
followed up on until very recently Il56ll57ll58ll . 

The work of Affleck et al. concerned the imaginary part of the scalar QED effective 
Lagrangian in a constant electric field. This quantity has been of much interest ever 
since Schwinger, building on earlier work by Sauter, Heisenberg, Euler, and Weisskopf 
[[59l [60l EEL showed that its existence implies the possibility of electron-positron pair 
creation in vacuum by the electric field ll62l . For small production rates this rate is simply 
given by twice the imaginary part itself, /production ~ 2Imj£f [2?]. Moreover, Schwinger 
was able to explicitly calculate the imaginary part in terms of a sum of exponentials, 



ImJ^ spinor (£) = ^3^ 2 L^exp 



Tin 
J 



(18) 



(j8 = eE/m 2 ). In this sum the fzth term relates to the coherent production of n pairs by 
the field. The appearance of E in the denominator of the exponents indicates that the 
pair creation effect is nonperturbative in nature. It also suggests an interpretation as a 
tunnel effect where a virtual electron-positron pair separates out along the field lines and 
extracts a sufficient amount of energy from the field to turn real. 

The corresponding formula for scalar QED differs only by a global factor and signs 
(due to the difference in statistics), 



Im^scalar(£) 



—J E— exp 



The production rates are exponentially small for 



Tin 

J 



(19) 



m 



£<£ crit = — =1.3xl0 18 V/m 
e 



(20) 



Until a few years ago producing an electric field close to this critical field strength E CY i t , 
and with a sufficient spatial extension, appeared far out of the reach of laboratory exper- 
iments. However, due to recent advances in laser technology it seems now conceivable 
that pair production could be seen in laser fields in the near future. The optical laser PO- 
LARIS, under construction at the Jena high-intensity laser facility, is projected to reach 
a maximal field strength of E mSLX ^2x 10 14 V/m [63], while the European X-ray free 
electron laser (XFEL), under construction at DESY, is expected to come even closer, 
E^ x » 1.2 x 10 16 V/m|l6l. 

For realistic laser experiments it is usually far from clear whether the constant field 
approximation is justified; it would be preferable to have generalizations of Schwinger's 



formula (18) to inhomogeneous and time-dependent fields. This is a subject which has 



been pursued by many authors, and many results have been obtained over the years (see, 
e.g. i3[MIE3[MI^[7^ However, considering this large volume of 

work there is a surprising lack of variety in the calculation methods used; apart from 
a few special field configurations, virtually all results of a more general nature have 
been obtained by WKB, or some variant of it. The formalism developed below is, while 
similar in spirit to WKB, technically quite different, and apparently more general. 

Let us start with retracing Affleck et al.'s lfT8ll recalculation of the Schwinger formula 
for scalar QED, ( [T9| ). At one loop, Feynman's representation ([!]) reads 



Tscalar [A] — 

Rescaling t = Tu, this becomes 



J~*L g _„? t J ^ xt -f T dz(i+ieA.x) (21) 



[A] = f~fz-™ 2T J®xA 



Tscalar 

The T integral has a stationary point at 



Y Jq dux 2 -\-ie Jq duA-x 



(22) 



yj f dux 2 



m 



(23) 



This allows us to calculate its imaginary part using a stationary phase approximation, 
yielding 



Imr sca i ar — 
with a new worldline action, 



1 Im [ 2)x \ I — e 
m J V T c 



27T — (myj J dux 2 -\-ie Jq 1 duA-x^j 



(24) 



S = m\l J dux 2 + ie J duA-x 



(25) 



We would like to calculate the remaining path integral using a stationary phase approxi- 



mation, too. The action (25) is stationary if 



m 



yj f dux 1 



ieF^yXy 



(26) 



Contracting with yields i 2 = constant = a 1 , and thus 



— ^2 



mx^ = ieaF^yXy 



(27) 



Thus the extremal action trajectory x cl (u), the "worldline instanton", will be a periodic 
solution of the Lorentz force equation, with a parameter a to be determined by the peri- 
odicity condition. Once the instanton is found, its worldline action immediately provides 
a semiclassical approximation for the imaginary part of the effective Lagrangian. Closer 
inspection shows that this approximation also corresponds to a weak field approximation 

ImJ^ scalar (£) E ~° e-^ cl l (28) 

For a the case of a constant electric field, E = (0, 0, E) = const., it is easy to see that the 
periodicity condition is solved by 

a n = — 2nn, n G Z + (29) 
eE 

The worldline instantons are simply circles in the t — z plane, with a winding number n: 



x n ( u ) — ~^ ^i,^2?cos(2n^M),sin(2n^w)J 
S[4] = nn — 

(30) 

The evaluation of the worldline action on the nth instanton thus yields just the nth ex- 
ponent in Schwinger's formula ( [T9| ). Moreover, Affleck et al. were able to compute also 
the prefactor, which here involves the determinant of fluctuations around the instanton 
path (see below). Thus, this approach provides a very simple and elegant rederivation of 
Schwinger's result for scalar QED. 

We will now sketch how to generalize this approach to general electric fields, as 
well as to the spinor QED case; see (56l |57l |58) for the details. Starting all over from 
Feynman's representation of the one-loop effective action in scalar QED ([T]), for the 
general case we prefer not to eliminate the T integral, but rather to first seek a stationary 
phase approximation for the path integral. The stationarity condition is again the Lorentz 
force equation, 



= 2ieF^ v (x)x v (31) 

We fix a point on the loop: 

[ ^x(T)e~ s ^= [d 4 x^ [ ^x(T)e~ s ^ (32) 

Jx(T)=x(0) J Jx(T)=x(0)=x(°) 

Let us assume that we have found a worldline instanton x c1 (t), a classical solution with 

x(T) = jc(0) = x(°\ We expand around x cl : 



Obtain the operator of quadratic fluctuations (Jacobi or Hessian matrix) A^y 



(33) 



where 



I d 2 d d 



(34) 



(X) 

Find the zero modes r/ v ; (t) of A^y, 



d 2 L 



dx^dxy ' 



R 



d 2 L 
dx^dxy 



(35) 



A M vT? v A) =0 



with initial value conditions 



(36) 



tj v A) (0) = 0, W(0) = 8 v?i , (ju,v = 1,2,3,4) 



(37) 



Evaluate them at T = T. Then, the final result for the semiclassical approximation 
becomes ll57ll 



Wt)^- 5 W t ) 
x(r)=x(o)=^°) 



ie P -s[i cl ](r) 



e e 



(4ttt) 2 



N 



det 



det 



(38) 



Note that the calculation of this quite general fluctuation determinant has been reduced to 
finding the determinant of the 4 x 4 matrix of zero modes. This remarkable simplification 
relies on a theorem by Levit and Smilansky [17511 . Finally, at the very end one does / dT 
using the stationary phase method again. 

For a number of classical special cases of "planar" non-constant fields, such as 
the single-pulse time dependent field [[68l|69l, the single-bump space dependent field 
Il67lh the sinusoidal time-dependence (66l[69l, and the sinusoidal space-dependence, the 
worldline instantons can be found explicitly in terms of special functions Il56lh leading 
to simple explicit formulas for the semiclassical exponent. For example, for the single- 
pulse field 



E{t) = £sech 2 (on) 



(39) 



the stationary action is ll56l 



m 2 7i ( 2 \ 
W = n—-[- 7== (40) 

where n is the winding number and y = ^ the "adiabaticity parameter" Il65ll . For the 
analogous single-bump field 



E(x 3 ) = £sech 2 (fot 3 ) (41) 

one finds 

m 2 7i ( 2 \ 
^bump = n— — (42) 

where y—mk/eE. Note that S pu i se decreases with 7, so that the pair production rate 
increases, while it is the other way round for Sbump- We believe that this is just an instance 
of a quite general fact: Inhomogeneity in time tends to shrink the size of the wordline 
instantons, leading to an increase in the pair production rate; spatial inhomogeneity 
increases the instanton size and decreases the pair production rate. 

The single-bump field also provides an excellent opportunity to test the validity of 
the semiclassical approximation, since for this case an exact integral representation has 
been obtained by Nikishov ll67lh suitable for numerical evaluation, and moreover there 
is a worldline Monte Carlo result ll50ll . As shown in figure |6j all three results are in close 
agreement over the whole range of the inhomogeneity parameter 7. 



Nikishov 

worldline instanton 

O Gies & Klingmiiller 



FIGURE 6. Plot of the imaginary part of the effective action for the field E(x) = £sech 2 (&x) as a 
function of the inhomogeneity parameter 7, normalized by the weak field limit of the "locally constant 
field approximation" (571 . The dotted line shows the result obtained from the worldline instanton approx- 
imation. The dashed line is the same ratio using a numerical integration of Nikishov 's exact expression 
1 67 ]. The circles represent the worldline Monte Carlo results of l50lL evaluated for eE/m 2 = 1. 



Note that the imaginary part vanishes for y > 1, which in the instanton approach 
simply means that instanton solutions cease to exist. Although mathematically this 
absence of instanton solutions allows one to conclude the vanishing of the imaginary part 
only in the semiclassical approximation, this example of the single-bump case supports 
the conjecture that it might, in fact, signal the complete absence of pair production. 

Although the method works very well for these classical cases, these are rather special 
configurations, and known to be amenable also to WKB methods. The existence of 
closed-form expressions for the worldline instanton can be expected only for a very 
restricted classes of fields. Much more interesting is the fact that the worldline instanton 
equations ( [27] ) and the zero mode equations ( [36] ) are ordinary differential equations, 
which leads us to expect that they can be solved numerically for more general classes of 
fields than have been treated by WKB. Dunne and Wang ll58ll have very recently applied 
this numerical approach to a class of electric fields which depend nontrivially on two 
spatial coordinates, parametrized by 



A 4 (x) 



(43) 



(y = mk/eE). Fig. 7^ shows their results for two examples, 



m 
m 



k(xi +X2) 
1+k^xj+xj) 

k(x l +x 2 )e- k2{x2 i +x ^ 



(44) 



O-O f(x,y)=(x+y)/(l+x +y )) 
& -H f(x,y)=(x+y)e x " y 




FIGURE 7. Pair creation rates for two nonplanar cases. 

Note that again Imr vanishes for large inhomogeneities, i.e. for fields of insufficient 
extent. 



Reproduced from l58l with the permission of the authors. 



Finally, making the transition from scalar to spinor QED in the instanton approach is 
straightforward if we use Feynman's original implementation of spin in the worldline 
path integral. Up to the global normalization, it amounts to multiplying with the spin 
factor ([4]) evaluated on the instanton trajectory, S[x cl ,A]. For the classical cases discussed 
above, and more generally for all "planar" fields the path ordering in the spin factor has 
no effect. Surprisingly, one obtains simply [f56l 



S[x c \A] 



4 cos 



eT / duE(x(u)) 



4(-ir 



(45) 



with n the winding number. Thus for planar fields one finds the same simple relation 
between IvnJ^^x^E) and Im^f sp i nor (£') as in the constant E case, eqs.([T8|),(19)! It is 
an interesting open question whether this property might even extend to generalelectric 
fields. 



5. SUMMARY 



The purpose of this short review was to show that the worldline approach in its various 
versions is turning into an efficient alternative to second-quantized methods for an 
increasing range of problems in QED. Unfortunately, despite of the restriction to QED 
it was not possible here to give due credit to all relevant work. Among other things, it 
was not possible here to discuss the worldline variational approach of J76l[73, nor the 
nonperturbative propagator calculations of ll78ll . 

Let me conclude with summarizing the main advantages which one can hope to 
achieve using the worldline formalism: (i) Compact parameter integral representations 
for arbitrary QED multiloop amplitudes (ii) Easy implementation of constant external 
fields (iii) Reliable numerical (Monte Carlo) results for one loop effective actions and 
Casimir energies (iv) Reliable numerical results for pair creation rates in arbitrary 
electric fields. 



REFERENCES 



1. R. P. Feynman, Phys. Rev. 80, 440 (1950) 

2. A. Morgan, Phys. Lett. B 351, 249 (1995), hep-p h/95Q2230| 

3. A. A. Migdal, Nucl. Phys. B 265, 594 (1986). 

4. C. D. Fosco, J. Sanchez-Guillen, and R. A.Vazquez, Phys. Rev. D 69, 105022 (2004), hep-th/03 10191 ; 
Phys. Rev. D 73, 045010 (2006), |hep-th/0511026| 

5. R. P. Feynman, Phys. Rev. 84, 108 (1951). 

6. E. S. Fradkin, Nucl. Phys. B 76, 588 (1966). 

7. F. A. Berezin, and M. S. Marinov, Ann. Phys. 104, 336 (1977). 

8. L. Brink, S. Deser, B. Zumino, P. Di Vecchia, and P. S. Howe, Phys. Lett. B 64, 435 (1976). 

9. A. Barducci, R. Casalbuoni, and L. Lusanna, Nuovo Cim. A 35, 377 (1976). 

10. L. Brink, P. Di Vecchia, and P. Howe, Nucl. Phys. B 118, 76 (1977). 



11. M.G. Schmidt, and C. Schubert, Phys. Rev. D 53, 2150 (1996), hep-th/9410100 

12. A. Strominger, Phys. Lett. B 101, 271 (1981). 

13. A. M. Polyakov, Mod. Phys. Lett. A 3, 335 (1988). 



14. H. Gies, and J. Hammerling, Phys. Rev. D 72, 065018 (2005), hep-th/0505072 

15. C. Schubert, Phys. Rep. 355, 73 (2001), [h^th/0101036| 



16. M. B. Halpern, and W. Siegel, Phys. Rev. D 16, 2486 (1977). 

17. M.B. Halpern, A. Jevicki, and P. Senjanoivc, Phys. Rev. D 16, 2476 (1977). 

18. I.K. Affleck, O. Alvarez, and N.S. Manton, Nucl. Phys. B 197, 509 (1982). 

19. Z. Bern, and D.A. Kosower, Nucl. Phys. B 379, 451 (1992). 

20. M.J. Strassler, Nucl. Phys. B 385, 145 (1992). 

21. M.G. Schmidt, and C. Schubert, Phys. Lett. B 331, 69 (1994), hep-th/ 9403158| 

22. K. Roland, and H.-T. S ato, Nucl. Phys. B 480, 99 (1996), [h^th/9604152| Nucl. Phys. B 515, 488 
(1998), hq>th/9709019| 

23. K. Daikouji, M. Shino, and Y. Sumino, Phys. Rev. D 53 , 4598 (1 996), |hep-ph/9508377| 

24. R. Shaisultanov, Phys. Lett. B 378, 354 (1996), hep-th/95 12142 



25. M. Reuter, M.G. Schmidt, and C. Schubert, Ann. Phys. (NY.) 259, 313 (1997), hep-th/9610191 

26. A. M. Polyakov, Gauge fields and strings, Harwood Academic Publishers, 1987. 

27. C. Schubert, Eur. Phys. J. C 5, 693 (1998), hep-th/97 10067 

28. P. Dai, and W. Siegel, YITP-SB-06-29, |hep-th/0608062| 

29. C. Itzykson, and J. Zuber, Quantum field theory, McGraw-Hill, 1985. 

30. D.J. Broadhurst, R. Delbourgo, and D. Kreimer, Phys. Lett. B 366, 421 (1996), piep^ph795 09296 

31. D.G.C. McKeon, and T.N. Sherry, Mod. Phys. Lett. A 9, 2167 (1994). 

32. S.L. Adler, and C. Schubert, Phys. Rev. Lett. 77, 1695 (1996), hep-th/9605035 

33. C. Schubert, Nucl. Phys. B 585, 407 (2000), hep-ph/0001288 



34. D. Fliegner, M. Reuter, M.G. Schmidt, and C. Schubert, Theor. Math. Phys. 113, 1442 (1997), |hep- 
lth/97041941 

35. B. Kors, and M.G. Schmidt, Eur. Phys, J. C 6, 175 (1999), |hepTh/9803144| 

36. G.V. Dunne, A. Huet, D. Rivera, and C. Schubert, JHEP 11, 013 (2006), hep-th/0609065 

37. G.V. Dunne, and C. Schubert, JHEP 0208, 0053 (2002), hep-th/0205004 

38. G.V. Dunne, and C. Schubert, JHEP 0206, 0042 (2002), hep-th/0205005j 

39. M.G. Schmidt, and C. Schubert, Phys. Lett. B 318, 438 (1993), hep-th/9309055 

40. D. Cangemi, E. D'Hoker, and G. Dunne, Phys. Rev. D 51, 2513 (1995), hep-th/9409113 

41. VP. Gusynin, and LA. Shovkovy, Can. J. Phys. 74, 282 (1996), hep-ph/9509383; /. Math. Phys. 
40, 5406 (1999), hep-th/9804143; LA. Shovkovy, AMS/IP Stud. Adv. Math. 13, 467 (1999), |hep] 
th/9902019, 

42. D.G.C. McKeon, and A. Rebhan, Phys. Rev. D 47, 5487 (1993), [hepTh/9211076| Phys. Rev. D 49, 
1047 (1994), hep-th/9306148] 

43. H.-T. Sato, /. Math. Phys. 40, 6407 (1999), hep-th/9809053| 

44. LA. Shovkovy, Phys. Lett. B 441, 313 (1998), hep-th/9806156 

45. D.G.C. McKeon, and A. Rebhan, Phys. Rev. D 48, 2891 (1993). 

46. A.I. Karanikas, and C.N . Ktorides, JHEP 9911, 033 (1999), hep-th/9905027 ; Phys. Lett. B 500, 75 
(2001),[hep-th/0008078l 

47. H. Gies, a nd K. Lang feld, Nucl. Phys. B 613, 353 (2001), [hep^ph /0 102 185; Int. J. Mod. Phys. A 17, 
966 (2002),|hep-ph/01 12198] 

48. M.G. Schmidt, and I.O. Stamatescu, hep-lat/0201002; Mod. Phys. Lett. A 18, 1499 (2003). 

49. H. Gies, K. Langfeld, and L. Moyaerts, JHEP 0306, 018 (2003), hep-th/0303264 

50. H. Gies, and K. Klingmuller, Phys. Rev. D 72, 065001 (2005), hep-ph/0505099 

51. H. Gies, J. Sanchez-Guillen, and R.A. Vazquez, JHEP 0508, 067 (200 5), | hep-th/0505275| 

52. H. Gies, and K. Klingmuller, /. Phys. A 39, 6415 (2006), [hl^tn/05 110921 

53. H. Gies, and K. Klingmuller, Phys. Rev. Lett. 96, 220401 (2006), quant-ph/0601094 

54. F. Bastianelli, O. Corradini, and P.A.G. Pisani, JHEP 0702, 059 (2007), hep-th/06 12236 

55. A. Bulgac, P. Magierski, and A. Wirzba, Phys. Rev. D 73, 025007 (2006), hep-th/05 11056 

56. G.V. Dunne, and C. Schubert, Phys. Rev. D 12, 105004 (2005), hep-th/0507174 



57. G.V. Dunne, H. Gies, C.Schubert, and Q.-h. Wang, Phys. Rev. D 73, 065028 (2006), hep-th/0602176 

58. G.V. Dunne, and Q.-h. Wang, hep-th/0608020 , 

59. F. Sauter, Zeit.f. Phys. 69, 742 (1931). 

60. W. Heisenberg, and H. Euler, Z. Phys. 98, 714 (1936). 

61. V. Weisskopf, Kong. Dans. Vid. Selsk. Math-fys. Medd. XIV 6 (1936). 

62. J. Schwinger, Phys. Rev. 82, 664 (1951). 

63. T. Heinzl et al., hep-ph/0601076 

64. TESLA technical design report. 



65. L.V. Keldysh, JETP 20, 1307 (1965). 

66. E. Brezin, and C. Itzykson, Phys. Rev. D 2, 1191 (1970). 

67. A.I. Nikishov, Nucl. Phys. B 21, 346 (1970). 

68. N.B. Narozhnyi, and A.I. Nikishov, Yad. Fiz. 11, 1072 (1970). 

69. V.S. Popov, JETP 34, 709 (1971); V.S. Popov, and M.S. Marinov, Yad. Fiz. 16, 809 (1972). 

70. M. Stone, Phys. Rev. D 14, 3568 (1976). 

71. W. Greiner, B. Miiller, and J. Rafelski, Quantum Electrodynamics of Strong Fields, Springer, 1985. 

72. A.B. Balantekin, J.E. Seger, and S.H. Fricke, BMP A 6, 695 (19912 

73. G.V. Dunne, and T. Hall, Phys. Rev. D 58, 105022 (1998 ),|hep-th/9807031| 

74. S.P. Kim, and D. Page, Phys. Rev. D 73, 065020 (2006), hep-th/0 30TT32| |hep-th/070T047 

75. S. Levit, and U. Smilansky, Ann. Phys. 103, 198 (1977). 

76. R. Rosenfelder, and A.W. Schreiber, Phys. Rev. D 53, 3337 (1996), [nucTtn7 9504002| Eur. Phys. J. C 
37, 161 (2004), hep-th/0406062] 

77. A. Alexandrou, R. Rosenfelder, and A.W. Schreiber, Phys. Rev. D 62, 085009 (2000), hep-th/0003253 

78. C. Savkli, J. Tjon, and F. Gross, Phys. Rev. C 60, 055210 (1999), hep-ph/9906211 ; Erratum-ibid. 61, 
069901 (2000); Phys. Rev. C 62, 116006 (2000),[hep-ph/99074451 



